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The structure of thin-film water on a BaF 2 (lll) surface under ambient conditions was studied using x-ray 
absorption spectroscopy from ambient to supercooled temperatures at relative humidity up to 95%. No 
hexagonal ice-like structure was observed in spite of the expected templating effect of the lattice-matched 
(111) surface. The oxygen K-edge x-ray absorption spectrum of liquid thin-film water on BaF 2 exhibits, at all 
temperatures, a strong resemblance to that of high-density phases for which the observed spectroscopic 
features correlate linearly with the density. Surprisingly, the highly compressed, high-density thin-film 
liquid water is found to be stable from ambient (300 K) to supercooled (259 K) temperatures, although a 
lower-density liquid would be expected at supercooled conditions. Molecular dynamics simulations indicate 
that the first layer water on BaF 2 ( 1 1 1 ) is indeed in a unique local structure that resembles high-density water, 
with a strongly collapsed second coordination shell. 



Water is unique in that it exists in several different identifiable amorphous forms, low-density amorph- 
ous (LDA), high-density amorphous (HDA), and very- high- density amorphous (VHDA), with a 
density difference of 33% between LDA and VHDA 1 " 4 . Several metastable amorphous forms have 
also been observed at various pressures and temperatures 5 ' 6 . It has been proposed that LDA and HDA are glassy 
forms of low-density liquid (LDL) and high-density liquid (HDL), respectively, in a liquid-liquid transition 
scenario 1 ' 7 " 10 , if these liquids would exist as distinct phases. There are many crystalline phases that resemble 
some of the amorphous forms, but much higher density ices can also be formed under pressure; e.g., ice VII has a 
density 74% higher than hexagonal ice I h u . All of these forms of crystalline ice phases formed at modestly high 
pressures, as well as amorphous or glassy water, typically have structures in which there is an inward collapse of 
the 2 nd coordination shell. Similarly, water under pressure at room temperature responds through a continuous 
collapse of the 2 nd shell with increasing pressure 12 where the extrapolated pair- correlation function (PCF) of HDL 
water has a broad 2 nd -shell distribution in between those of HDA and VHDA 2 ; a similar structural rearrangement 
has been observed in simulations of water in silica nano-confmement 13 . Water thus seems to have the ability to 
exhibit rather dramatic structural differences with a large range of densities and only minor differences in 
enthalpy 1 . Here we demonstrate that it is possible to generate an extremely compressed HDL-like form of water 
at ambient conditions at the BaF 2 (lll)/water interface, which acts as the pressurizing medium. 

The influence of the substrate on interfacial water has often been discussed in terms of ice-like and liquid-like 
structural signatures 14 " 17 . In this respect, BaF 2 (lll) terminated by F"-Ba 2+ -F" trilayers, with a surface lattice 
parameter of 4.39 ±0.1 A, is very close to that of the basal plane of ice I h (4.52 A). Thus, an epitaxial layer of ice 
growth is expected to form on BaF 2 (l 1 1) at ambient conditions. Although earlier optical absorption experiments 
suggested the formation of ice-like overlayers under low temperature 18 and ambient conditions 19 , theoretical 
estimates of the stability of the epitaxial ice layers are not consistent with these experimental observations 20 . Also 
in nano-confmement Raman scattering experiments on water have found a non-freezing water layer at the water/ 
solid interface down to extremely low temperatures 21 . Arguments about the stability of epitaxial wetting layers 
based on molecular dynamics simulations 22 have not been corroborated because of the effects of surface 
defects 23 " 25 . Theoretical calculations 20 using the ASP-W2K potential find an adsorption energy of —39.78 kj/ 
mol for an isolated water molecule on BaF 2 . Higher coverages resulted in a bilayer hexagonal ice film with an 
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increase of the adsorption energy to —51.4 kj/mol, which, albeit 
higher than the enthalpy of condensation of bulk water ( — 44 kj/ 
mol), demonstrates only a weak trend towards wetting. 

In the present study, the molecular structure of thin-film water 
under ambient conditions (up to 2.4 monolayers (ML)) as well as ice 
under low-temperature conditions in ultra-high vacuum (UHV) 
were elucidated on BaF 2 (lll) using x-ray absorption spectroscopy 
(XAS) and x-ray photoelectron spectroscopy (XPS). O K-edge XAS is 
a powerful local probe for identifying the structures of water in solid 
and liquid phases; the spectral features are directly related to the 
hydrogen-bonding environments 26 . The experimental data are com- 
plemented with molecular dynamics (MD) simulations, where the 
first interfacial layer indicates a very high local density (1.23 g/cm 3 ) 
and the O-O pair- correlation function shows a near- complete col- 
lapse of the 2 nd shell, well beyond what is seen in HDL and close to 
merging with the first shell to form a two-dimensional hexagonal 
structure. 

Results 

In order to characterize the structure of the thin water films we used 
O K-edge XAS to obtain spectral signatures indicative of specific 
hydrogen-bonding configurations. The obtained spectra of thin 
layers of liquid water, however, do not exhibit the expected templat- 
ing effect resulting in hexagonal ice, but instead show strong sim- 
ilarity to high-pressure ices, such as ice VII with a density of 1.6 g/ 
cm 3 . XAS spectra measured at 1.5 Torr water at 30% (288 K), 60% 
(270 K), and 90% (259 K) relative humidity (RH), with the polar- 
ization vector of the incident x-ray out-of-plane and in -plane, are 
shown in Figs. l(a and 1(b)), respectively. These films correspond 
to 1.2, 1.8, and 2.4 ML coverage of water (see Supplementary 
Information for details). At all RHs studied, the XAS spectral features 
are rather similar: strong pre-edge and weak post-edge features are 
observed, indicating that, irrespective of the polarization vector ori- 
entation and the temperature (—300-259 K), the structure of thin- 
film water on BaF 2 (l 1 1) does not depend on coverage. At 30% RH, a 
monolayer of water forms, as seen in the uptake measurements. 

The pre-edge feature at 535 eV has been proposed to be mainly 
related to weakened H -bonds 26,27 whereas the post- edge feature 
appearing at 540-542 eV is due to H-bonds and is particularly strong 
for tetrahedrally H-bonded water as in hexagonal ice (ref. 26,28 and 




530 535 540 545 

Photon energy (eV) 



530 535 540 545 550 

Photon energy (eV) 



Figure 1 | XAS spectra of water on BaF 2 (lll) at 1.5 Torr water pressure 
with the polarization vector (a) perpendicular and (b) parallel to the 
surface plane. RH was varied from 30% to 90%. The incident angle (a) of 
in-plane polarized light and the take-off angle (6) are —10° and 35°, 
respectively. 



discussions therein regarding various interpretations). The main- 
edge feature, at —537 eV is enhanced upon formation of HDA 29 
and various crystalline high-pressure ice phases such as I n , Iyi> Ivn> 
and I V iii 30 - This enhancement is related to the fact that the second 
shell collapses from the tetrahedral distance of 4.52 A in ice I h to 
much shorter distances, leading to an overall higher density; note 
that this also slightly increases the pre-edge intensity. 

In Fig. 2a, we compare the XAS spectrum of thin-film water on 
BaF 2 (lll) at 90% RH with the spectrum of crystalline ice formed at 
low pressures (Ih or I c ) on BaF 2 (lll) and with the relative change 
between ice I h and I V n 30 . The latter corresponds to a high-pressure 
crystalline ice phase with a density of 1.60 g/cm 3 . First, we note that 
the thin water film XAS spectrum lacks the post-edge feature which is 
seen in crystalline ice. Furthermore no low pressure ice-like spectral 
features were observed at any coverage as RH was increased from 
30% to 90%, even at temperatures as low as 259 K. Considering the 
arguments of surface-templated ice formation, the templating effect 
of BaF 2 (lll) on the structure of the water film is expected to be 
strongest under this condition. Secondly, we instead note that the 
water thin-film spectrum has a close resemblance with the high- 
pressure spectrum of ice I V n- 

In Fig. 2b, we note that the XAS spectrum of thin-film water on 
BaF 2 (lll) at ambient conditions more closely resembles the spec- 
trum of liquid water than that of low-density crystalline ice on BaF 2 , 
as seen from the weak post-edge and strong main-edge and pre-edge 
features. Because the post-edge is enhanced for tetrahedral H-bond- 
ing configurations, we infer that low- density crystalline ice-like local 
arrangements are very sparse in the water layer on BaF 2 (lll). This 
inference is quite surprising because the substrate temperature of 
259 K is deep in the supercooled state of water, and the driving force 
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Figure 2 | (a) Comparison of XAS spectra of (i) bulk ice I h (dashed), ice VII 
(I V ii) (solid) 30 with the spectra of (ii) 2 ML crystalline ice and (iii) 
supercooled water on BaF 2 (lll) at 1.5 Torr at 259 K (90% RH). 
Crystalline ice on BaF 2 (lll) was prepared by annealing 10 ML LDA to 
175 K. XAS spectra of crystalline ice on BaF 2 (lll) were recorded via 
exciting with in-plane polarized light, (b) Comparison of XAS spectra of 
(i) supercooled water on BaF 2 (lll) at 1.5 Torr at 259 K (90% RH) with 
(ii, solid) liquid water, and (ii, dashed) 6 M NaCl solution measured at 
room temperature 31 . The supercooled water spectrum is an average of two 
spectra obtained using in-plane and out-of-plane polarizations, shown in 
Fig. 1. The shaded area represents the difference between bulk ambient 
liquid and thin-film water. XAS spectra of liquid water and the 6 M NaCl 
solution were recorded in transmission mode 31 . The sample was embedded 
between two Si 3 N 4 windows separated by a 300 nm thick spacer. All 
spectra are normalized by area. 



SCIENTIFIC REPORTS | 3 : 1074 | DOI: 1 0.1 038/srep01 074 



2 



for tetrahedral structuring should be rather strong according to 
expectations from thermodynamics arguments that enthalpy and 
entropy of water films converge to values typical of ice 19 . Contrary 
to these enthalpy and entropy estimates from quantification of IR 
spectral features, the tetrahedral structure of ice for water films on 
BaF 2 (lll) has not been found to be the most favorable based on 
theoretical calculations 20 and polarizability measurements 25 . 

Even more surprising is the observation that, due to the loss of 
intensity in the post-edge region, the spectral width spanned by the 
main-edge and post-edge peaks of supercooled water on BaF 2 (l 1 1) is 
somewhat smaller than that for liquid water. In Fig. 2b we also 
compare the thin-film water spectrum with the spectrum of 6 M 
NaCl solution at room temperature 31 and note in the small post-edge 
a stronger similarity for thin-film liquid on BaF 2 (l 1 1) with the NaCl 




Figure 3 | Oxygen-oxygen pair-correlation functions (a) calculated 
within the first 1 A layer of H 2 0 on the BaF 2 (lll) surface from MD 
simulations, full line T=260 K, dotted line T=290 K; (b) of HDL 34 ; (c) of 
a NaCl solution of concentration 1:10 salt to water mole ratio 44 ; 
(d) calculated within a 1 A layer of bulk H 2 0 from MD simulation at 
T=260 K; (e) calculated within the second 1 A layer of H 2 0 on the 
BaF 2 (lll) surface from MD simulations, full line T=260 K, dotted line 
T=290 K. Note that the pair-correlation functions calculated within a 1 A 
layer of water (a, d, e) are normalized differently than in (b, c), cf. 
Supplementary Information. Inset: Histogram of distances between the 
oxygen atoms of H 2 0 and the surface atoms. 



solution than with bulk water. In alkaline aqueous solutions the 
interaction of water with the alkali cations is weak 26,31 , while neutron 
diffraction measurements 32 indicate that the resulting structure is 
still similar to that of water under pressure, which is also supported 
by MD simulations 33 . The close spectral similarity of the thin-film 
water spectrum and the spectra of the NaCl solution and I V n indi- 
cates that the water structure should resemble more a HDL-related 
structure as typically seen upon applying pressure on water with an 
inward collapse of the 2 nd shell. 

From the MD simulations, we find a strong layering of water 
molecules parallel to the surface. The inset of Fig. 3 shows a dramatic 
localization of the first water layer; the distribution of the molecules 
within subsequent layers is broader, but the influence of the 1 st layer 
on the 2 nd and 3 rd layers is evident. To further determine the structure 
within these layers, we calculated oxygen-oxygen PCFs within layers 
of defined thickness perpendicular to the surface and compared the 
calculated PCFs with the experimental PCFs of HDL water 34 , NaCl 
solution 31 and the PCF of a slab of bulk TIP4P/2005 water of the same 
thickness as for water on BaF 2 (Fig. 3). The first three peaks in the 
TIP4P/2005 bulk water PCF represent the 1 st , 2 nd , and 3 rd coordina- 
tion shells located at 2.75, 4.4, and 6.7 A, respectively, in this model. 
The 1 st layer of 1 A thickness of thin-film TIP4P/2005 water on 
BaF 2 (l 1 1), on the other hand, shows a significantly altered structure 
compared to the 2 nd and 3 rd layers of the same thickness. The 2 nd 
coordination shell of the 1 st layer is collapsed as in the case of 
HDL 34 , and the 3 rd peak shifted inwards to 6 A, indicating that the 
3 rd coordination shell adapts to the high- density profile. The 2 nd and 
3 rd layers of the thin film, on the other hand, show a structure similar 
to that of the bulk. For T=260 K, the structure is slightly more 
pronounced than for T=290 K, but does not change qualitatively. 
These findings are universal and do not change with the different 
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Figure 4 | A snapshot of the first water layer on BaF 2 (lll). F" and Ba 2+ 

ions are color-coded by purple and green, respectively. Red and white 
sticks represent oxygen and hydrogen atoms of water molecules. 
Intermolecular hydrogen bonding and water-BaF 2 (lll) bonding are 
indicated by dashed lines. H-bond definition: O-O distance < 3 A, O-O-H 
angle <30 deg. 
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force fields commonly utilized to simulate the properties of water 
(see Supplementary Information for details). 

In Fig. 4 we show a snapshot from the simulation with a strong 
predominance of chain-like conformations in the layer in contact 
with the substrate. In comparison to low density crystalline ice, the 
structure of the 1 st layer is more similar to close packing at the 
expense of breaking H-bonds. The density in the 1 st layer is estimated 
to be 1.23 g/cm 3 (see Supplementary Information for details), which 
is similar to the density of VHD A; this similarity also fits well with the 
2 nd -shell distance of 3.3 A 1 ' 4 . However, there are more severe distor- 
tions from tetrahedral coordination in the 1 st shell in comparison to 
VHDA which we attribute to the thermal fluctuations prominent 
near room temperature. 

Discussion 

Previous surface x-ray scattering and diffraction studies of the water- 
silica 35 and water-Cu(lll) 36 interfaces have shown that the first 
water layer has a structure with higher density than that of either 
ice or liquid water. A high-density form of interfacial water has also 
been identified for water on NaCl(OOl) in MD simulations 37 as well as 
for water in nano- confinement 13 ' 21 . In a recent study of the role of 
lattice mismatch in ice nucleation Cox et a/. 38 found that the angular 
flexibility in the second shell allowed closer packing at their model 
interface. These findings can now be interpreted in the context of the 
present results with the formation of a highly compressed, HDL-like 
phase at the interface. In fact, the HDL-like phase persists well below 
freezing without any signatures of ice formation in a way similar to 
the freezing point lowering seen in NaCl solutions 39 . 

We have shown that the BaF 2 (lll) surface does not act as an ice 
formation template for ambient and supercooled water films, not- 
withstanding the close lattice match. The XAS spectrum of water on 
BaF 2 (lll), on the contrary, indicates limited tetrahedral H-bonding 
at 300-259 K, where the similarity to spectra of high density ices 
suggests a highly compressed HDL-like water structure even in the 
supercooled regime. The O-O pair-correlation function of the first 
water layer from MD simulations fully supports this high-density 
form of water at the interface. 

Methods 

Measurement techniques. Relative humidity- dependent water uptake curves were 
obtained using XPS to monitor the water coverage (see Supplementary Information 
for details). The BaF 2 (lll) surface was exposed to 1.5 Torr water vapor and the 
uptake was carefully controlled by reducing the sample temperature in an 
incremental fashion. 2.4 ML of water was the maximum coverage at 90% RH and 
uptake curves are in good agreement with the ones measured using IR spectroscopy 19 . 

The low temperature, UHV experiments on ice on BaF 2 (lll) were performed at 
the soft x-ray Spectroscopy Beamline (BL 13.2) at the Stanford Synchrotron Radiation 
Lightsource (SSRL). X-ray absorption spectra were acquired by the total electron yield 
(TEY) method. Ambient pressure XPS and XAS experiments were performed at the 
Molecular Environmental Science Beamline (BL 1 1.0.2) at the Advanced Light Source 
(ALS). The differentially pumped electron energy spectrometer (Specs Phobios 150) 
enables measurements at pressures up to 5 Torr 40 . XAS data at 1.5 Torr H 2 0 were 
acquired by the Auger electron yield (AEY) method using a kinetic energy detection 
window of 370-400 eV. Details of the sample preparation and data analysis are 
presented in figure captions and in the Supplementary Information. 

MD simulations. MD simulations were performed to model the behavior of H 2 0 on a 
fluorine-terminated BaF 2 (l 1 1) surface (256 BaF 2 , 500 H 2 0). We used the GROMACS 
package ver. 4.5.3 41 , the OPLS-AA force field 42 , and the TIP4P/2005 water model 43 ; 
results using alternative force-fields are discussed in the Supplementary Information. 
Each simulation was equilibrated for at least 5 ns, generic thermodynamic properties 
were monitored to ensure equilibration. The production runs were performed in the 
NPT ensemble for 20 ns at a time step of 1 fs using velocity rescaling to set the 
temperature; the Berendsen barostat was applied for the pressure. All simulation 
results shown here were obtained at 1 bar and 290 K or 260 K, respectively. 
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